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\^ , We investigate the pattern of asymmetries in the colour-octet model 

I for inclusive production of charmonium with polarised initial state 

• hadrons and photons. We find that singly polarised asymmetries van- 

5^. ish. The double polarised asymmetry for J/ip involves a new combi- 

QhI nation of long-distance matrix element. We also investigate the asym- 

metries for r]c, Xc aud T production. The asymmetries distinguish 
well between different models of charmonium production. 
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Abstract 



1 Introduction 



A large volume of data on the inclusive hadro-production of J / ip has a good 
phenomenological description in terms of the non-relativistic QCD (NRQCD) 
based colour octet model |l|, ^ |^ . There is no proof of NRQCD factorisation 
for these cross sections and the agreement between data and model may be 
considered phenomenological. Since other production mechanisms for this 
cross section have also been postulated ^ ^ ^, it is useful to explore the 
pattern of spin asymmetries in this model. Polarisation asymmetries for 
inclusive charmonium production in various models have been investigated 
earlier 0, ^, 0. Here we perform a detailed analysis in the context of the 
colour-octet model. In the near future, experiments at HERA and RHIC can 
use these predictions to assess the relative importance of colour octet and 
other contributions to inclusive J/ip production. 

The physics of heavy quarks can be expected to be dominated by its 
(large) mass, m, and hence computable in perturbative QCD. However, for 
quarkonium systems the scales mv (relative momentum) and mv"^ (energy) 
are also important, since the quarks are expected to be non-relativistic in 
the rest frame of the quarkonium. An organised way to include the effects of 
the dimensionless parameter v (relative velocity) is provided by an effective 
low-energy field theory called non-relativistic QCD (NRQCD) A proof 
of factorisation can be given for several processes |1TT| . 

In this framework, a heavy quark pair created by a short- distance process 
binds into a quarkonium state over scales that are longer by powers of 1/v. 
Assuming factorisation, the inclusive cross-section for the production of a 
meson H can be written as 

MH) = \dT Y. ^JyOiC'^'LA) (1) 

n={a,S,L,J) '"' 

where $ is a flux factor, dT is the phase space measure and the C's are short 
distance coefficients for the production of a heavy quark pair of given colour 
(a), spin (S) and orbital angular momentum (L) coupled to a total angular 
momentum (J). The naive scaling dimension of the matrix element is dn- It 
describes the long-distance physics involved in the transition of the pair to 
the observed hadron H. The colour octet model gets its name from the fact 
that the colour index can run over singlet as well as octet representations. 
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The C's have the usual perturbation expansion in powers of the strong 
couphng as- The matrix elements {O^ (^^~^^Lj)) scale as y3+2L+2_E+4M -^j^Qj-e 
E and M is the minimum number of gluon electric and magnetic dipole 
transitions connecting the produced pair to the target hadron. Consequently, 
the sum in eq. (|l]) is a double power series expansion in v and a^. In this 
paper we consider all terms of order a^f " (n < 7) for hadro-production and 
order aagv"' {n < 7) for (almost elastic) photo-production. Lepto-production 
will be considered elsewhere. 

Most of the matrix elements {O^ i^^^^Lj)) required for the spin-averaged 
J/ip cross section to this order are now being extracted from data; some 
of them are known with fair accuracy. They were first obtained from the 
phenomenology of P-state charmonium production at large Pt in Tevatron 
data [1^, where they are required for consistency of the perturbation series 
[ P^ . The matrix elements for S-state hadrons were also seen to be important 
for phenomenology A global analysis of Tevatron data has now been 



performed and many of these matrix elements have been fixed A 
linear combination of octet matrix elements which appears for the first time in 
inclusive J /ip production has been estimated through an analysis of (almost) 
elastic photo-production data p!7| . 

In section 2 we write down the coefficient functions appearing in eq. (|^) 
for the production of various resonances with polarised initial states. The 
spin asymmetries are written down in section 3. This section also contains 
a discussion of the kind of the phenomenological signals expected from the 
colour octet model and experimental checks on competing models. The con- 
cluding section contains a discussion of the effects of higher orders in v. 



2 The Coefficient Functions 

In the past, the coefficient functions were computed using a covariant projec- 
tion technique 0, |18|, A new algorithm for computing these coefficient 



functions from the hard-scattering diagrams has recently been developed 
We follow this technique but use spherical tensor methods for the simplifica- 
tion of angular momentum projections. Our phase conventions are given in 
Appendix 0. 

The heavy quark {Q) and anti-quark [Q) have momenta p and p respec- 
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tively (p^ = = m?). These are written as 

p'^ = ^P'^ + Lfg* and p^" = ^P" - ^q', (2) 

where P is the centre of mass (CM) momentum of the pair and q is the relative 
momentum in the CM frame. This has only spatial (Euclidean) components 
(denoted by Latin indices) and is boosted by the Lorentz transformation 
with mixed Lorentz (in Greek) and spatial indices. The latter contract over 
spatial indices of any space-like object and boosts it properly, provided it has 
vanishing time component in the CM frame of the pair. 

In this paper we are concerned only with 2 2 processes giving a QQ 
pair in the final state. If the initial 4-momenta are taken to be ki and k2, 
then P = ki + k2- Further, K {= ki — /C2) is spacelike and orthogonal to 
P and hence is also boosted properly by . Finally, for convenience, we 
choose the direction of K (in the rest frame of the pair) to be the z direction 
and quantise the angular momenta along this direction. In a hadron-hadron 
collision this corresponds to choosing the common direction of ki, ^2 and P 
to be the z direction. 

The perturbative QCD cross sections are written down as usual and then 
matched to eq. (|l|) to extract the coefficient functions. In order to do this, 
the Dirac spinors for Q and Q are written in terms of Pauli spinors {tp for 
Q and x for Q), and the matrix element is Taylor expanded in q. Since we 
only consider terms upto order f ^, we can drop higher than linear terms in 
q. The angular momentum projections are all effected by recouplings of the 
spin-1 spherical tensor . The matrix elements for the processes gg QQ 
and qq QQ have been constructed and used to compute the coefficient 
functions for the cross sections with unpolarised initial states. These provide 
checks on our construction. 

2.1 qq QQ 

Keeping all terms upto linear order in g, we find that the colour averaged 
squared matrix element for the lowest order qq process is given by 

|^|2 ^ 477^^^ _ ^^^^ _ .^.^^ {^ia^Ty^^a^T'^x), (3) 

where h and h' are the helicities of the two quarks, a denotes a Pauli matrix 
and T a colour generator. The terms linear in h and h' couple to the operator 
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ieijk{x^ cr^ipip'^ cr^ x) ■ Since this is a negative parity operator, it does not appear 
in the NRQCD action. As a result the single spin asymmetry vanishes. We 
discuss this point later. 

After inserting the proper projection operator, ^.^{M), for the final state 
hadron H with helicity M, and converting Euclidean indices to helicity in- 
dices (A and A'), we obtain the expression 



(4) 



A summation over A and A' is understood. When the helicity of the final 
state H is observed, well-known angular momentum techniques can be used 
to reduce the matrix element to its standard form given in eq. (|lD. For 
the J/ip, the results of are reproduced. If the helicity is unobserved, 
then a sum over M can be performed, and the matrix element reduces to 

(-i)X-A'(0re^i))/3. 

The last ingredients in the calculation of the coefficient functions are the 
phase space measure and the flux factor — 

where s = P"^ = {ki + ^2)^. The coefficient functions for the unpolarised 
cross section (C) and double polarised asymmetries (C) can now be written 
down. The results are — 

ci's!) = -CCS!) = (6) 



C has been obtained previously in [^. Note that these coefficient functions 
are to be multiplied by the dimension 3 operator matrix element {Og {^Si)) . 
The parton level double polarised asymmetry is —1 independent of the final 
hadron. 



2.2 gg^QQ 

We turn now to the gg QQ processes. In a ghost-free gauge there are three 
of these — an s-channel gluon exchange diagram and the t and u channel 
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heavy quark exchange diagrams. The summed matrix element breaks up 
into three colour structures 



A 



g'e';,{h,h)et{h,h') 



-Sabs'"" + -d^'^'D^" + -f^b^F^^" 
6 2 2 . 



(7) 



Note that Yang's theorem prevents the amplitudes S and D from coupling to 
a spin-1 state No such restriction applies to the amplitude F. Note also, 
that this theorem restricts only the state of the QQ pair. The transition from 
this to the final observed hadron may involve further soft gluon radiation. 
As a result, the angular momentum of the final hadron is not restricted by 
Yang's theorem. 

We take the density matrix for any initial state gluon of momentum k 
and helicity h to be given by 



e^'{k]h)e*''{k]h) 



1 
5 



-9 



k^'p" + k'^p'' p^k^k"" 



k-P 



(k-py 



ih 
■2FrP 



(8) 



This gauge choice avoids ghosts due to the three gluon vertex. 

The square of the amplitude S after colour averaging is given by 



1^1 



2 2 

36rr? 



(1 + hh')i^{x^^UH{M)ij^x) + ^ {h,2Sx',-2 + 6x,-2Sy,2) 



m 



X^[aq]lijUH{M)i:^[aq]l,x 



1-hh') (x^[aqfo^UH{M)^naq]^oX 



(9) 

We have introduced the notation [abl^ to denote two 3- vectors, a and b, 
coupled to a spherical tensor with angular momentum J and helicity A. In 
this form the result is remarkably simple, compared to the same expressions 
written out in Euclidean components. The square of the D amplitude is the 
same, apart from colour generators inserted between the Pauli spinors. The 
square of the amphtude F is simply 



o 2 2 

Sir 



' ' 32m" 
with a helicity index on q. 



f{l + hh') {x^q'^TyUHiM)^^q^T''x 



(10) 
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In eq. (^) we have dropped terms linear in h and h\ which would con- 
tribute to the single spin asymmetries. These contain a non-zero coefficient 
multiplying the operator 

xVn(i/)V^V* [- '-D^^ x) + he. (11) 

Since this is a negative parity operator, it does not appear in the NRQCD 
action. As a result, single spin asymmetries vanish. The same argument 
shows the absence of single spin asymmetries in |-Dp. 

This mechanism is quite general. The single spin asymmetries are pro- 
jected by squaring the perturbative amplitude and contracting one set of 
indices with the symmetric part of eq. (H) and another with the antisym- 
metric part. This would always couple to a negative parity operator. Such 
operators cannot appear in the NRQCD action. Consequently single spin 
asymmetries would vanish in all processes involving only strong or electro- 
magnetic vertices. Such an argument can also be developed for qq initial 
states. 

The L = (S-state) and L = 1 (P-state) matrix elements have mass 
dimensions 3 and 5 respectively. The coefficient functions can then be written 
down — 

ci'sl) = CCS',) = 



9m 

2 2 

vr a 



3m 



s 



ccsl) = ccsl) = fc'i'si), (12) 

CCPl) = -C^PI) = fc^^Pi), 

Im^ 

All the parton level asymmetries in gg —>■ QQ are exactly 1 except the 
^P2 asymmetry, which is —1. This fact has interesting phenomenological 
consequences. 



ccpf) = ccpf) = 
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2.3 -fg QQ 

The coefficient functions can be easily written down for photo-production of 
J/ip in the almost elastic limit z ^ 1. The inelasticity parameter z is defined 
by 

(13) 



p ■ 

where p, k^ and k^ are the 4-momenta of the target proton, final state char- 
monium and the photon respectively. The only amplitude appearing in this 
process is proportional to the amplitude D defined in eq. (|^). To linear or- 
der in q there are three coefficient functions — C{^Sq), C{^Pq) and C(^P|). 
These differ from the corresponding coefficients in eq. ( [T^ ) by an overall 
common factor of 64a/15as- Note that these formulae can only be applied 
after subtracting the diffractive component at z = 1 p2| and summing the 
remainder over a small range of z (say z > 0.95). 



3 Asymmetries 

In this section we consider the asymmetries for J/ip and other charmonium 
states produced in polarised hadron-hadron or photon-hadron cross sections. 
The quantities 

^ = zEh,h'^ih,h'), 

(14) 

Aa = \j:h,h'hh'aih,h'), 

computed to order f ^ are listed in Appendix Polarisation asymmetries 
are given by 

Parton level cross sections and asymmetries are distinguished from the phys- 
ical (hadronic) quantities by a caret over any symbol for the former. 

The hadronic cross sections and asymmetries involve parton distributions 
in the projectile (P) and target (T) hadrons. These are usually combined 
into a parton luminosity factor which depends on whether one studies the 
integrated cross section or distributions in the Feynman variable Xp. For 
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unpolarised initial states, the luminosities relevant to Xp distributions are — 



Cg — I x^gp(x^) x^Qj-^x^), 

(16) 



xl + At 



Y.f x+qfp{x+)x-qf,{x-) + {P ^ T) 



where the sum in the definition of Cq is over fiavours. Here r = 4m^/ S, where 
S is the CM energy of the colliding hadrons. These luminosities multiply the 
cross sections in eq. (^). We have used the notation 

1 



x± = -{^xI + At±Xp). (17) 

The luminosities appropriate for total cross sections are given by integrals 
over Xp. For polarised scatterings, the parton distributions are replaced 
by the polarised parton distributions Ag{x) and Aq{x). These polarised 
luminosities are denoted by ACg and ACg. 

The hadronic cross sections in photo-production of charmonium involve 
the polarised or unpolarised gluon densities at 

where m^v is the nucleon mass and is the photon energy in the rest frame 
of the target nucleon. 

3.1 "yp ^ J/i/j 

The polarisation asymmetry for photo-produced J/ ip is given by 

where the combination of matrix elements G and G are given in eq. (|37|). 
Since G is unknown, it is not possible to predict this asymmetry. Photo- 
production experiments must measure G. Present knowledge (summarised 
in Appendix 0) indicates that A['yp — ^ J/ '4^) could be larger than the gluon 
asymmetry Agjg^ and have the same sign as the latter. 
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Upto order f xj states are not produced in 7p reactions, hence the 
asymmetries A{'-)p xj) ^''^^ undefined at this order. For the same reason, 
AIjp — > rjc) is undefined at order v^. 



3.2 pp ^ J/l/j 

The colour octet model predictions for the polarisation asymmetries for xo, 
Xi and r]c are 











y4(pp - 


^ Xo) 



(20) 



We can use the relations in eq. (|^) to simplify the results. Since {Oi'' i^Pj)) 
0.32±0.04 GeV^ and {Of{^Si)) = 0.0098±0.0013 GeV^ [|T5|, and £g > £ 



the qq processes can be neglected to a good approximation. As a result we 
find 

MpP Xo) ~ -A{pp X2) ~ ^^g/^g = A{pp r/J. (21) 

The asymmetry for J / ip production can also be written down. Again, the 
matrix elements for the qq subprocesses are much smaller than those for the 
gg subprocesses. As a result, to a good approximation. 



A{pp — > J l^j) 
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5m^9 + 15Po (Of (^Pq) ) - AB, {Of (^P^) ) ™, 

5m^e + 15Po(Of e^)) + 4P2(Of ('P2)) ^ ^' 

(22) 

where Pj is the branching ratio for the decay of xj to J/ip. If can be 
measured in 7p reactions, then this asymmetry is predicted without free 
parameters. 



It should be understood that in eqs. (19-22), the ratio of matrix elements 
is a kinematics-independent factor. Thus, the same ratio appears for inclu- 
sive integrated asymmetries and for x^-dependence of the asymmetries. The 
kinematics is subsumed into the parton density factors. As a result, apart 
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from the cross checks between hadro- and photo-production, each of these 
has internal cross checks in the agreement between inclusive asymmetries and 
their x^-dependence. 



3.3 T production 

It is interesting to note the differences in T production asymmetries. Since 
f ^ is much smaller for the bottomonium system than for charmonium, and 
the branching fractions for Xb — *■ 7^" are larger, the factors corresponding to 
and B in the analogue of eq. (p2D are relatively small. Consequently, 



The colour octet components in this asymmetry are small, and the dominant 
part is the same as in the colour singlet model. 

In 'jp —>■ T reactions, however, the colour octet parts (of order v"^) are 
dominant, since there are no colour singlet parts at lower order in v. Thus, 
the prediction for the polarisation asymmetry is 



Ahp^m = ^ I (24) 




where the subscript b on the matrix elements denote that they are taken for 
the bottomonium system. 

3.4 Distinction between models 

It is known that the colour octet contributions to the integrated J ftp cross 
section is important for unpolarised scattering. The normalisation of the 



cross section requires the matrix element defined in eq. (pTf). If these 
contributions were absent, as in the colour- singlet model Q], then, since Bq 
is small, 

AC 

A{pp —>■ J/ip) ^ A{pp —>■ X2) = -pr^ (colour singlet model). (25) 

Falsification of this relation by the data would indicate the necessity of colour 
octet components in the cross section. Detailed predictions based on the 
colour singlet model have been extensively explored 
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An alternative model for the production of charmonium is the semi-local 
duality model [^]. In this model the cross section for producing any charmo- 
nium resonance is obtained by taking the sum of the qq QQ and gg QQ 
cross sections (integrated in the range 4m^ < 5" < 4M|,, where M^) is the 
mass of the D meson) and multiplying them by a phenomenological number 
derived from the total cross sections. This model is difficult to distinguish 
from the colour-singlet model using data on unpolarised cross sections. How- 
ever, the semi- local duality model predicts 

A{^p ^ H) 
A{pp H) 

independent of the final state hadron H [§]. This is very different from the 
prediction of the colour octet model, and an experimental test should be 
simple. 

Higher twist effects ^ give rise to more complicated Xp dependence of 
cross sections. Large effects of this kind would violate the Xp dependences 
imphcit in eqs. (|1)-(|2|). 

4 Discussion 

In section 3 we have considered all subprocesses upto order w^. However, 
from the analysis of data on the unpolarised cross sections, it is known that 
certain aspects of the phenomenology require matrix elements upto order f ^ 
[0]. These will be considered in detail elsewhere. However, it is possible to 
identify the major qualitative changes when these effects are included. 

The most important change is in the production of xj- Upto order 
these resonances are produced only through the qq subprocess and the am- 
plitude S in the gg subprocess. As a result, upto this order, A{^p xj) is 
undefined. At order w^, xj can be produced through the D amplitude in the 
gg subprocess, for example by the matrix element {O^'' Sq)) . Intermediate 
pair states with L = 2 also occur. Consequently all these states can also 
be produced in 7p collisions. Therefore, A{'jp — > xj) becomes non-zero and 
gets related to Ag/g. 



Agjx) 

^^"^^ . (semi — local duality), (26) 

r -L r — 
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At order t>'^, A{pp xo, X2) is related to ACg/Cg, whereas A{pp xi) = 
ACg/Cq. At order v^, because of the processes discussed above, A{pp — ^ xi) 
is well approximated by ACg/Cg. 

These changes in the xj cross sections and asymmetries reflect on those 
for the J/ip. Final state J/ip's can be produced by radiative decays of these 
additional Xj's- Eqs. (|T9|) and ( p2D will then involve new combinations of 
matrix elements. 

Some of these considerations are very well illustrated by rjc production at 
order v'^. At the leading order v^, A{pp — > rjc) is given by the expression in 
eq. (pOl). There is no rjc production in 7p collisions at this order. At order 
v'^, however, 

Hl9 - Vc) = S{s - ^m')^^^ {ore So)). (27) 
As a result, we obtain the asymmetry 

^h,^,J = M£). (28) 

9{x) 

In pp collisions, order f ^ terms are obtained through the D amplitude which 
gives the same matrix element as in eq. p7| ) above. In addition, the ma- 
trix element (Og^l^Pi)) enters through the F amplitude in gg — » QQ. The 
asymmetry becomes 

A( ^ ^ 3ACg^-64Cg{OfCS,)) 

^ '''^ ^ACg^ + MC,{ai^('S{)) ' (29) 

where vl> = 8(Of (^^o)) + lh{Ot{^So)) + ^{Ofi^P,)). 

Since \E' is of order whereas {O^i^So)) is of order t>^, one expects about 10% 
corrections to the leading order result (eq. ^). This is a marginal effect, but 
may be observable. Note that the linear combination ^ involves the matrix 
element {(Dr{}Pi)), from the amplitude F, not seen in 77^ production in 7p 
collisions. 

For the bottomonium system, the substantially smaller value of f ^ implies 
that these higher order corrections in v are rather small. 

The remaining part of the analysis presented in this paper is robust; the 
results do not change as one takes higher orders in v into account. We 
conclude by summarising the most robust predictions of the colour octet 
model. 
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All single spin asymmetries vanish in processes which do not involve at 
least one W oi Z. 



Polarisation asymmetries for inclusive production and the Xp distribu- 
tion of any charmonium state involve the same combination of non- 
perturbative matrix element. The difference is only in the parton dis- 
tribution factors. 

Double spin asymmetries for xo, X2 and rjc are simply related to the 
gluon asymmetries in hadro-production (eq. ^1]). The asymmetry for 
X2 is (almost) equal and opposite to those for xo and rjc. 

The J/ ip asymmetries are simply related to gluon asymmetries in both 
hadro- and photo-production (eqs. |19|, 



The charmonium asymmetries provide easy distinction between differ- 



ent models of charmonium production (eqs. ^) 



The asymmetry in the hadro-production of T is very similiar to the 
results of the colour singlet model (eq. ^3|). That for photo-produced 
T is similiar to the J/ip (eq. 



A Phase Conventions for Spherical Tensors 

Conversion of Euclidean indices to helicity labels is effected by a unitary 
transformation. This can be specified by the spherical components of a gen- 
eral vector V — 

V, = -^{V^ + tVy), v., = -^{-V,, + tVy), Vo = V,. (30) 

With this phase convention, we may define the transformation of the Eu- 
clidean metric — 

5,, ^ (-1)"5a,-a'. (31) 

Since the metric is not positive definite, it makes sense to define covariant 
and contravariant indices. We choose the components of a covariant vec- 
tor to satisfy eq. (^), and hence transform as bras under rotations. The 
contravariant vector then transforms as a ket. 
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Clebsch-Gordan coefficients define unitary transformations between bases 
of quantum states specified by two (commuting) angular momenta ji and j2 
and tlieir resultant J, and is simply the matrix element (ji, Tni;j2, m2\J, M). 
Standard orthogonality relations follow from the completeness of either ba- 
sis. Following pSf we choose these matrix elements to be real and fix 



(ii) Ji! J25 j2bi+j25 ji+j2) = 1- Spherical tensors are irreducible components 
of Euclidean tensors under rotations and can be constructed by ffist convert- 
ing Euclidean indices to helicity labels and then recoupling them through the 
Clebsch-Gordan coefficients. 

The metric can be extended to higher tensors, and is easily written in 
terms of a Clebsch-Gordan coefficient — 



(-1)'5a,-a' = ^{2J+l){J\-J\'\m). (32) 

The Levi-Civita tensor can also be expressed in terms of a Clebsch-Gordan 
coefficient — 

e^Jk (-l)^(l,A;l,/i'|l,-z/). (33) 

These two tensors suffice for the problem at hand. 

Simple angular momentum recouplings give the leading order heavy-quark 
spin symmetry results 

{OU'Pj)) = (2J+l)(Of,8(^Po)), (34) 
{On^S,)) = (2J+l)(Of(3^i)). 

Order f ^ corrections to these are harder to compute. 



B Cross Sections 

In this appendix we list the parton level cross sections and polarisation asym- 
metries for charmonium production to order for rjc, order for xj cind 
order for J/ip. The results for unpolarised scattering have been obtained 
before pl 

For the J/ip, the ^Sq and the ^Pj matrix elements are of order v"^ . In 
addition, the ^Sf matrix element from the qq process also contributes to the 
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same order. All other matrix elements are of higher order. However, the J /il) 
final state can also be obtained by radiative decays of the xj states. These 
arise at order through the ^P] matrix elements in the gg subprocess and 
at the same order through the matrix element from the qq subprocess. 
Upto order these are the only channels which have to be summed to 
obtain the cross section for J /if) in hadro-production. For the rjc the 
matrix element in the gg subprocess is the leading term. Other terms are 
subleading at least by v^. 

H, 13— 



The unpolarised cross sections are well-known |I8| 



^{99 W) 



©."♦(■So 



c^i.99 - 


-^Vc) = 


6{s 


^{99 - 


^ Xo) = 


6{s 


(^{99 - 


^ X2) = 


S{s 


a{qq - 


-H) = 


6{s 



Am 



m 
s 2 

2^ TT ag 



36m^ 

„3„ 2 



Pi 



(35) 

In the qq process the hadron H can be taken to be anything. The cross 
section for direct J/?/' production involves the linear combination of matrix 
elements determined in 0, . 

The polarisation asymmetries are given by Aa/a. The numerators are 
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listed here — 



Aa{gg - 




Aa{gg - 


^ Xo) 


Aa{gg - 


^ X2) 


Aa{qq - 





6{s — Aw? 
5{s — 4m^ 
— — Aw? 
—5{s — Aw 



36rr? 
3 
w 



+ 



5m 



P2 



36m' 



12m' 
3 



|(Of(3Po)) 



3n,2 
S- 



27 w^ 



r)H(3 

^8 



Si] 



(36) 

The polarisation asymmetry for direct J/ ip production contains a new hnear 
combination of matrix elements, as can be seen from the expression above. 
The expressions for 'jg —>■ J/ip cross sections can be obtained by multiplying 
those for gg — >• J/ip by 6Aa/15as. 

Heavy quark spin symmetry can be used to reduce the matrix elements in 
polarised and unpolarised direct J/?/' production to the two combinations — 



e 
e 



So] 



7 

w 

1 

w 



(unpolarised), 
(polarised). 



(37) 



There have been attempts to extract the combination G from experiments; 
it is found that the value hes in the range 0.02 and 0.03 [0, |1^. Another 
combination. 



+ 



m 



Pn 



(38) 



has also been extracted, and its value lies between 0.04 and 0.066 |T3|, ITB]. 
The uncertainties reflect the spread in data as well as the choice of parton 
densities, and perhaps also indicate the need for higher order computations. 
If these two combinations are used to determine = 2^9 — B, then we find 
that G/G lies between 1.7 and 6.6. This spread is larger than the order 
errors involved in writing eq. (pTj). In a previous work it has been assumed, 
following a suggestion in [jl6[, that {O^^^ {^Pq)) ^ 0, and hence G/G f» 1. 
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